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■ We solve a problem, which while not fitting into the usual paradigm, can be viewed 

■ as a quantum computation. Suppose we are given a quantum system described by an 
N dimensional Hilbert space with a Hamiltonian of the form E\w){w\ where \w) is an 
unknown (normalized) state. We show how to discover \w) by adding a Hamiltonian 
(independent of \w)) and evolving for a time proportional to N 1 / 2 /E. We show that 
this time is optimally short. This process is an analog analogue to Grover's algorithm, 
a computation on a conventional (!) quantum computer which locates a marked item 
from an unsorted list of N items in a number of steps proportional to iV 1 / 2 . 
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Although a quantum computer, beyond certain elementary gates, has not yet been con- 
structed, a paradigm |L] for quantum computation is in place. A quantum computer is 
envisaged as acting on a collection of spin 1/2 particles sitting at specified sites. Each ele- 
mentary operation is a unitary transformation which acts on the spins at one or two sites. 
A quantum computer program, or algorithm, is a definite sequence of such unitary transfor- 
mations. For a given initial spin state, the output of the program is the spin state after the 
sequence of transformations has acted. The length of the algorithm is equal to the number 
of elementary unitary transformations which make up the algorithm. 

This framework for quantum computation is general enough that any ordinary digital 
computer program can be turned into a quantum computer algorithm. (It is required that 
the ordinary program be reversible; however any ordinary computer program can be written 
in reversible code.) Quantum computers can go beyond ordinary computers when they act on 
superpositions of states and take advantage of interference effects. An example of a quantum 
algorithm which outperforms any classical algorithm designed to solve the same problem is 
the Grover algorithm 0. Here we are given a function f(a) defined on the integers a from 1 
to N. The function has the property that it takes the value 1 on just a single element of its 
domain, w, and it has the value for all a ^ w. With only the ability to call the function 
/, the task is to find w. On a classical computer this requires, on average, N/2 calls of the 
function /. However Grover showed that with a quantum computer w can be found with 
of order N 1 / 2 function calls. This remarkable speed-up illustrates the power of quantum 
computation. (In the appendix we explain how the Grover algorithm works.) 

In this paper we consider quantum computation differently, as controlled Hamiltonian 
time evolution of a system, obeying the Schrodinger equation 

ij t m=H(tm, a) 

which is designed to solve a specified problem. We illustrate this with an example. Suppose 
we are given a Hamiltonian in an N dimensional vector space and we are told that the 
Hamiltonian has one eigenvalue E ^ and all the others are 0. The task is to find the 
eigenvector \w) which has eigenvalue E. We now give a solution to this problem and then 
explain in what sense it is optimal. 

We are given 

H w = E\w)(w\ (2) 

with \w) unspecified and (w\w) = 1. Pick some normalized vector \s) which of course does 
not depend on \w) since we don't yet know what \w) is. Now add to H w the "driving" 
Hamiltonian 

H D = E\s)(s\ (3) 

so that the full Hamiltonian is 

H = H W + H D . (4) 
We now calculate the time evolution of the state \ip w ,t) which at t — is \s), 

\^ W) t) = e~ mt \s). (5) 
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It suffices to confine our attention to the two dimensional subspace spanned by \s) and 
\w). The vectors \s) and \w) are (generally) not orthogonal and we call their inner product 

x, 

(s\w) = x (6) 

where x can be taken to be real and positive since any phase in (s\w) can ultimately be 
absorbed in \s). We will discuss the expected size of x shortly. Now the vectors 

1 



|r) 



(\s) — x\w)) 



and \w) are orthonormal. In the \w), \r) basis the Hamiltonian 



(7) 



is 



H = E 



1 +x 2 
x^fl 



x c 



X 



and 



X 



x- 



Now a simple calculation gives 



-iEt 



x cos(Ext) — i sin(Ext) 



\/\ — x 2 cos(Ext) 

Thus we see that at time t the probability of finding the state \w) is 

P(t) = sin 2 (Ext) + x 2 cos 2 (Ext) 
and that at a time t m given by 



7T 



2Ex 



(9) 
(10) 

(11) 
(12) 



the probability is one. 



How big do we expect x to be? In an iV dimensional complex vector space, if you pick 
two normalized vectors at random (uniformly on the 2N-1 dimensional unit sphere), then the 
expected value of the inner product squared is 1/N so we know that the expected value of x 
is of order N~ l l 2 . Thus starting with \s), for the probability of finding \w) to be appreciable 
we must wait a time of order N l l 2 / E. This is the analog analogue of the Grover algorithm 
result. 

Note that the eigenvalues of the Hamiltonian (§) are E(l ± x). Thus the difference in 
eigenvalues is (2xE) which is of order E/N 1 / 2 . By the time-energy uncertainty principle, 
the time required to evolve substantially, that is from \s) to \w), must be of order N l l 2 /E 
which is the time we found. You might think that by increasing the energy difference, that 
is for example, by using H D = E'\s)(s\ with E' ^> E you could speed up the procedure for 
finding \w). However the next result shows that this is not the case. 

We now show that our procedure for finding \w), in a time which grows like N 1 / 2 /E, is 
optimally short. The proof we give here is the analog analogue of the oracle proof [[| which 
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can be used to show that the Grover algorithm is optimal for the problem it sets out to solve. 
Again we are given the Hamiltonian H w = E\w)(w\ and we wish to add some Hamiltonian 
Ho(t) to it which drives the system to a state which allows us to determine \w). In an iV 
dimensional vector space, there are N linearly independent choices for \w). We can pick 
these to be a basis for the vector space and we then have 

]T H w = EY J \w)(w\ = E. (13) 

w w 



The idea of the proof is this: Start with some initial |u>)-independent state \i) and evolve 
it with the Hamiltonian 

H = H w + H D (t). (14) 

After a time t the state we get must be substantially different from what we would have 
gotten using H w , + Ho(t) or else we can not tell \w) from \w'). Let 



i-n\i>w,t) = (H w + H D (t))\ip w ,t) 
at 



(15) 



with 



\ip w ,0) = \i). 

In order for \ip w , t) to differ sufficiently from \ip w/ , t) it is certainly necessary that, for all (but 
one) w, \ip w ,t) differs sufficiently from any |«;)-independent vector. (If some of the \ip w ,t) 
were very close to a particular | ^-independent vector, we could not tell them apart.) Let 
\if),t) evolve with H D (t), that is, 



ij t \^,t) = H D (t)\ij),t) 



(16) 



with 



foM) = |i). 

We will use \ip,t) as a |w)-independent vector which the \ipw,t) must differ from. We require 

2 

t to be large enough that \ip w ,t) — \ip,t) > e for some fixed e which implies 



Y}\^ w ,t)-\^,t) 



> Ne. 



(17) 



Now consider 



\ip w ,t) - \iP,t)\\ 2 = -2Re j t ^j W} t\^,t) 



dt 



(18) 



which upon using ( P~5| ) and (|TBJ) gives 



dt 



\ip w ,t) - \if>,t) 



= 2 lm(il> w ,t\H w \il>,t) 
< 2\(^ W} t\H w \^ } t}\ 



(19) 



< 2 



H w \ip,t) 
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We now sum on w and use the fact that if Yli=i \ a i\ 2 — 1 then Y^Li |ai| < N 1 / 2 along with 
(|T3|) to obtain 

(f Elll^*)-^,*) 2 < 2 ^ 1/2 - (20) 



dt 

Since 0) = 0) we have 



£||V«»*>-hM>|| < 2EN 1 'H. (21) 

Therefore in order to satisfy (|17|) we must have 

eN 1 / 2 

«>- 51 - (22) 

This shows that the we have chosen allows us to determine \w) as quickly as possible in 
terms of N. 

Appendix: The Grover Algorithm 

We are given a function f(a) with a = 1,. . .N such that f(w) = 1 and f(a) = for 
a 7^ w. We assume that the function f(a) can be calculated using ordinary (reversible) 
computer code. The goal is to find w. Classically this requires, on average, N/2 evaluations 
of the function /. 

We now explain how the Grover algorithm solves this problem; see also [|J. The quantum 
computer acts on a vector space which has an orthonormal basis \a) with a — 1, ... N. It is 
possible to write a quantum computer algorithm which implements the unitary transforma- 
tion 

U f \a) = (-l) /(a) |a)- (Al) 

Equivalently we can write 

U f = l-2\w)(w\ (A2) 

The quantum computer algorithm which implements Uf requires two evaluations of the 
function / because it is necessary to erase certain work bits which we have supressed. It is 
also assumed that the ordinary code which is used to evaluate / has a length which does 
not grow like to a positive power. Then the number of two bit quantum computer steps 
required to evaluate / will also not grow as fast as iV to a power. 

Now consider the vector 

l«> = ^EI«>- (^) 

1 " a 

It is also possible to write quantum computer code which implements the unitary operator 

Z7» = 2|s><s| — 1. (A4) 

The number of two bit operations required to implement U s grows more slowly than to 
any positive power. 
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The Grover algorithm consists of letting the operator U s Uf act k times on the vector \s). 
To see what happens we can restrict our attention to the two dimensional subspace spanned 
by \s) and \w). Let 



so that \w) and |r) form an orthonormal basis for the relevant subspace. In the \w), \r) basis 
the operator U s Uf takes the form 



where cos 9 



U s U f = 

2/N. This implies that 



cos 9 
sin 9 



— sin 9 

cos 9 



cos{k9) — sin(k9) 
sin(k9) cos(k9) 



(A6) 



(A7) 



Now for N large 9 ~ 2N l / 2 so each application of U s Uf is a rotation by an angle ~ 2N 1 / 2 . 
In the \w), \r) basis, the initial state \s) is 

JV-V2 



\*) 



:i - ^) 1/2 



(as 



which is very close to |r). However after k steps where k9 = n/2 the algorithm has rotated 
the initial state to lie (almost) along \w). This requires k ~ nN 1 ^ 2 /4 steps. Each step 
actually requires two evaluations of / so the number of evaluations of / required to find w 
grows like N 1 ^ 2 . Accordingly the number of two bit operations required to implement the 
algorithm also grows like N 1 / 2 . 
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